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(54) Title: METHOD FOR ATTENUATION OF REVERBERATIONS USING A PRESSURE-VELOCITY BOTTOM CABLE 
(57) Abstract 



A method is described that attenuates reverberations in ocean bouom cable 
data used to generate a seismic profile. First, a seismic wave Is generated and 
recorded by pressure and particle velocity sensor pairs. Adding the product of 
the pressure data times the absolute value of the velocity data and the product 
of the velocity-data times the absolute value of the pressure data eliminates- 
reverberations, leaving a combined signal that is a function of only the primary 
wave response signal. Dividing this combined signal by 2 and a sensitivity seating 
factor, and then taking the square root of the result produces the magnitude of the 
primary signal. Replacing the positive/negative sign of the combined signal that 
existed before taking the absolute value provides the proper polarity. The data are 
then incorporated into a seismic profile. 
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METHOD FOR ATTENUATION OF REVERBERATIONS 
USING A PRESSURE-VELOCITY BOTTOM CABLE 

Field of the Invention 

This invention relates to seismic profiling of formations located under bodies of 
water using pressure-velocity bottom cables. 

Background of the Invention 

Marine seismic exploration is an extremely important tool for location of off- 
shore reserves. One known procedure for marine seismic exploration involves use of an 
ocean bottom cable as illustrated in prior art Fig. 1. Surveys using ocean bottom cables are 
typically employed in areas populated with numerous obstacles, such as drilling and 
production platforms. In this technique, several miles of bottom cables 1 1 (only one shown 
in the Fig. 1) are deployed along the sea floor 13 by vessel 15, Usually, multiple cables 11 
are deployed in parallel, as shown in Fig. 2. The bottom cable 11 is provided with a 
plurality of sensor pairs 17 placed at regular intervals along the cable, each sensor pair 17 
containing a pressure sensor (e.g., hydrophone) and particle velocity sensor (e.g., a 
geophone). Acoustic energy is generated in the vicinity of the cable using an air gun array 
or a marine vibrator array 19. The source wavelet travels downward through the earth and 
is partially reflected by subsurface layers (formation 21 in Fig. 1) that present an acoustic 
impedance contrast. The primary reflected wavelet 23 travels upwardly from the 
subsurface layer, and the pressure waves generated by the upward-traveling reflection are 
detected by the sensor pairs 17. 

Seismic exploration using ocean bottom cables is complicated by secondary 
waves such as wave 25, known as "ghosts," that are received by the sensors pairs 17 as 
downward-traveling reflections after reflecting off the air/ water boundary at the surface 29. 
The air/water boundary is an efficient reflector, and thus the ghosts are significant in 
amplitude and are difficult to differentiate from the primary waves. These ghosts adversely 
affect the data obtained during the exploration by attenuating certain frequencies. In 
addition, when the water depth is large, the spectral ghost notches fall in the seismic 
frequency band and drastically affect the seismic resolution. Resolution is further 
complicated by multiple reflection waves and water layer reverberations such as wave 27. 

The purpose of using both hydrophones and geophones in the ocean bottom 
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cable is to capitalize on the differences between these two types of sensors to attenuate the 
downgoing waves which include the ghosts and the water layer reverberations. Their 
responses to the primary reflections are in phase, but are 180" out of phase to the ghosts 
and to the reverberations. 

Fig. 3a illustrates amplitude versus time hydrophone response at the water 
bottom. For the hydrophone, at time ti (the primary wave), the hydrophone response is 
defined as T. At time t2 (the first water layer reverberation), the hydrophone response is 
-(l+r)T; at time t3 (the second reverberation) . it is r(l + r)T; and at time t4(the third 
reverberation) . it is -r2(|+r)T, where r is between 0 and 1 and is the water bottom 
reflectivity at each receiver position. Additional reverberations continue to decrease in 
amplitude. 

The geophone amplitude versus time response is shown in Fig. 3b. The 
amplitude at ti is MT. where M is a sensitivity scaling factor that depends upon the 
particular type of sensors used. At time t2, the geophone response is (l-r)MT. At time t3i 
the response is -r(l-r)MT, and at time 14 is r2(l-r)MT. 

It is apparent from Figs. 3a and 3b that while the primary geophone and 
hydrophone responses are in phase, the responses to the water layer reverberations are 180° 
out of phase. Thus, the attenuation of the reverberations can be achieved by adding the 
hydrophone and the geophone signals together after the signals have been suitably scaled. 
Theoretically, the scale factor S = (l+r)/(l-r), where r is the water bottom reflectivity, 
should be applied to the geophone data. Determination of the water bottom reflectivity 
coefficient r depends upon the acoustic impedance of the bottom material. Thus, the scale 
factor S can vary among different sensor pair locations on the same cable. 

There are several known methods for deriving the scaling factors for geophone 
signals. U.S. Patent No.5,235.554 describes a method where a calibration survey is used 
to estimate the water bottom reflection coefficient. In such a calibration survey, a low 
energy source is fired over each sensor pair, and the scale is determined from the ratio of 
the peaks of the first arrivals of the hydrophone and geophone signals. Collection of this 
survey data requires additional time and cost over and above the data acquisition phase of 
the survey. U.S. Patent No.5,396,472 describes a method to derive the water bottom 
reflection coefficient that eliminates the need for a separate calibration survey, but involves 
complex mathematics including summing the pressure and velocity signals, multiplying the 
results by the inverse Backus operator, and then solving for the water bottom reflectivity r 
using an optimization algorithm. U.S. Patent No. 5,365.492 describes a method wherein 
the hydrophone signal is used first to adaptively remove noise from the geophone signal. 
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and then the cleaned geophone signals are scaled by a scaling factor and added to the 
hydrophone signals. The resulting signal is then auto-correlated and the relative amplitude 
of the first auto-correlation's function-side lobe is measured. The optimum scale factor for 
the geophone data is then found by optimizing the value of the scale factor with respect to 
the first side-lobe amplitude of the auto-correlation. 

Summary of the Invention 

In one aspect, the invention relates to a method of generating a seismic profile 
using an ocean bottom cable having at least one sensor pair comprising a pressure sensor 
and a velocity sensor, comprising generating a seismic wave for reflection from earth 
strata; collecting pressure data with the pressure sensor and velocity data with the velocity 
sensor; multiplying the pressure data with the absolute value of the velocity data to produce 
a first result; multiplying the velocity data with the absolute value of the pressure data to 
produce a second result; summing the first result and the second result to create a third 
result; dividing the third result by a factor of 2 to obtain a fourth result; dividing the fourth 
result by a sensitivity scaling factor; recording a positive/negative sign for the fourth result; 
taking the square root of the absolute value of the fourth result to produce a fifth result; 
replacing the positive/negative sign of the fourth result into the fifth result to produce a 
sixth result; and incorporating the sixth result into the seismic profile. 

In another aspect, the invention relates to a method of generating a seismic 
profile using an ocean bottom cable having a first sensor pair and a second sensor pair, 
each of the first sensor pair and the second sensor pair comprising a pressure sensor and a 
velocity sensor, comprising generating a seismic wave for reflection from earth strata; 
collecting pressure data with the pressure sensors and velocity data with the velocity 
sensors; multiplying the pressure data from the first sensor pair with the absolute value of 
the velocity data from the first sensor pair to produce a first result; multiplying the velocity 
data from the first sensor pair with the absolute value of the pressure data from the first 
sensor pair to produce a second result; summing the first result and the second result to 
create a third result; dividing the third result by a factor of 2 to obtain a fourth result; 
dividing the fourth result by a sensitivity scaling factor; recording a positive/negative sign 
for the fourth result; taking the square root of the absolute value of the fourth result to 
produce a fifth result; replacing the positive/negative sign of the fourth result into the fifth 
result to produce a sixth result; multiplying the pressure data from the second sensor pair 
with the absolute value of the velocity data from the second sensor pair to produce a 
seventh result; multiplying the velocity data from the second sensor pair with the absolute 
value of the pressure data from the second sensor pair to produce an eighth result; summing 
the seventh result and the eighth result to create a ninth result; dividing the ninth result by a 
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factor of 2 to obtain a tenth result; dividing the tenth result by a sensitivity scaling factor; 
recording a positive/negative sign for the tenth result; taking the square root of the absolute 
value of the tenth result to produce a eleventh result; replacing the positive/negative sign of 
the tenth result into the eleventh result to produce a twelfth result; and incorporating the 
sixth result and the twelfth result into the seismic profile. 

In another aspect, the invention relates to a method of attenuating reverberations 
in seismic data, the seismic data including pressure daU and velocity data containing 
reverberations, comprising reading the pressure data from at least one pressure sensor; 
reading the velocity data from at least one velocity sensor; and processing the pressure data 
and the velocity data to produce composite data wherein the reverberations have been 
attenuated, comprising the steps of: multiplying the pressure data with the absolute value of 
the velocity data to produce a first result; multiplying the velocity data with the absolute 
value of the pressure data to produce a second result; summing the first result and the 
second result to create a third result; dividing the third result by a factor of 2 to obtain a 
fourth result; dividing the fourth result by a sensitivity scaling factor; recording a 
positive/negative sign for the fourth result; taking the square root of the absolute value of 
the fourth result to produce a fifth result; and replacing the positive/negative sign of the 
fourth result into the fifth result to produce said composite data. 

Brief Description of the Drawings 

Fig. 1 is an illustration of a conventional ocean bottom cable technique; 

Fig. 2 is an illustration of an array of ocean bottom cables parallel to each 

other; 

Figs. 3a and 3b are graphs illustrating hydrophone and geophpne response at the 
water bottom; 

Fig. 4 is a flow chart illustrating a method in accordance with an embodiment of 
the invention; 

Figs. 5a-5c are graphs illustrating results obtained from the responses shown in 
Figs. 3a and 3b using a method in accordance with the invention; 

Fig. 6 is a graph of hydrophone-geophone primary reflection data with 
reverberations; and 
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Fig. 7 is a graph of the data of Fig. 6 after application of a method in 
accordance with the invention. 

Description of the Preferred Embodiments 

. , Preferred embodiments of the invention will now be described with reference to 
the accompanying drawings. 

The invention allows unwanted reverberations in ocean bottom cable (OBC) data 
to be attenuated without the need for determining water bottom reflectivity and without 
requiring complex mathematical operations. Data may be collected by any known method 
including, e.g., the method shown in prior art Figs. 1 and 2. 

Fig. 4 is a flow chart of the steps employed to attenuate these reverberations in 
accordance with an embodiment of the invention. First, a seismic wave (51) is created and 
data are recorded with pressure and particle velocity sensors (53). Once the data are 
collected, the pressure response signal is multiplied by the absolute value of the velocity 
signal (55) at each time t, producing a result as shown in Fig. 5a. The combined primary 
wave response signal at time t] now has an amplitude of MT^; at t2, the amplitude of the 
signal is -M(l-r2)T2; at t3, the amplitude of the signal is Mr2(l-r2)T2; and finally, at t4. 
the amplitude is equal to -Mt^(\-t^)T^ (additional reverberations are not shown). 

Similarly, the response signal from the velocity sensor is multiplied by the 
absolute value of the pressure signal (57) at each value of time t. Fig. 5b shows the result 
of this operation. Combined primary wave response signal at time has an amplitude of 

Mt2, identical to the amplitude of the primary wave response signal in Fig. 5a. At time 
t2, the amplitude of the combined response signal for the first reverberation is equal to 

M(l-r2)T2; at 13, the amplitude of the combined response signal is now -Mr2(l-r2)T2; at 

t4, the amplitude is Mr4(Nr2)T2. 

As can be seen from Figs. 5a and 5b, the multiplication process has made all of 
the reverberations have identical amplitudes, but opposite polarities as shown by their 
positive/negative signs. Adding these signals at each time interval (59) produces a graph 
with only one response signal located at tj with a magnitude of 2MT^. This result is 

shown in Fig. 5c. This resultant signal is a function of the primary signal alone; all 
reverberations have been cancelled out by their respective counterparts. This signal is then 
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divided by two (61) and by the sensitivity scaling factor M (63). The sensitivity scaling 
factor M is dependent only on the type of sensors used in the bottom cable and this is 
specified by the manufacturer of the ocean bottom cable. Taking the square root of the 
absolute value of this last result (65) produces the true amplitude of the primary signal. 
Replacing the sign that existed before taking the absolute value in the previous step (67) 
provides the proper polarity for the signal. This result is then incorporated into a seismic 
profile (69). 

Fig. 6 illustrates synthetic OBC data containing six clusters of data (73, 75, 77, 
79, and 81). Each cluster includes one or more primary reflections 83 and a plurality of 
reverberations 85. It is clear from the data shown in Fig. 6 that the reverberations 85 are 
significant in amplitude and will seriously distort any seismic profile produced using this 
data if they are not successfully attenuated. 

Fig. 7 illustrates the data of Fig. 6 after application of the method described 
above and illustrated in Fig, 4. As is clearly shown, only the primary waves 83 arc 
remaining in the data of Fig. 7, Thus, much higher resolution can be achieved in the 
seismic profile, without necessitating expensive operations to determine a scaling factor for 
each receiver position in the survey. Once this data is produced with the reverberations 
attenuated, the processing of the seismic data is straightforward. The processing sequence 
by which such data as the data shown in Fig, 7 are converted in a final seismic profile that 
can be geologically interpretated is well known in the art and thus will not be described in 
detail herein. For instance, this can include deconvolution, static corrections, velocity 
determination, normal moveout corrections, dip moveout corrections, stack and migration. 
Detailed description of seismic data processing can be found in any of numerous texts on 
the subject, including for instance, "Seismic data processing" by Oz Yilmaz (Society of 
Exploration Geophysicists, Tulsa, Oklahoma, 1987). 

Various embodiments of the invention have been shown and described. 
However, the invention is not so limited, but rather is limited only the scope of the 
amended claims. 
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CLAIMS 

1. A method of generating a seismic profile using an ocean bottom cable having at 
least one sensor pair comprising a pressure sensor and a velocity sensor, comprising: 

generating a seismic wave for reflection from earth strata; 

collecting pressure data with the pressure sensor and velocity dala with the 
velocity sensor; 

multiplying the pressure data with the absolute value of the velocity data to 
produce a first result; 

multiplying the velocity data with the absolute value of the pressure data to 
produce a second result; 

summing the first result and the second result to create a third result; 

dividing the third result by a factor of 2 to obtain a fourth result; 

dividing the fourth result by a sensitivity scaling factor; 

recording a positive/negative sign for the fourth result; 

taking the square root of the absolute value of the fourth result to produce a fifth 

result; 

replacing the positive/negative sign of the fourth result into the fifth result to 
produce a sixth result; and 

incorporating the sixth result into the seismic profile. 

2. The method of claim 1, wherein the velocity sensor is a geophone. 

3. The method of claim 1 or claim 2, wherein the pressure sensor is a hydrophone. 

4. A method of generating a seismic profile using an ocean bottom cable having a 
first sensor pair and a second sensor pair, each of the first sensor pair and the second 
sensor pair comprising a pressure sensor and a velocity sensor, comprising: 

generating a seismic wave for refieciion from earth strata; 
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collecting pressure data with the pressure sensors and velocity data with the 
velocity sensors; 

multiplying the pressure data from the first sensor pair with the absolute value of 
the velocity data from the first sensor pair to produce a first result; 

multiplying the velocity data from the first sensor pair with the absolute value of 
the pressure data from the first sensor pair to produce a second result; 

summing the first result and the second result to create a third result; 

dividing the third result by a factor of 2 to obtain a fourth result; 

dividing the fourth result by a sensitivity scaling factor; 

recording a positive/negative sign for the fourth result; 

taking the square root of the absolute value of the fourth result to produce a fifth 

result; 

replacing the positive/negative sign of the fourth result into the fifth result to 
produce a sixth result; 

multiplying the pressure data from the second sensor pair with the absolute value 
of the velocity data from the second sensor pair to produce a seventh result; 

multiplying the velocity data from the second sensor pair with the absolute value 
of the pressure data from the second sensor pair to produce an eighth result; 

summing the seventh result and the eighth result to create a ninth result; 

dividing the ninth result by a factor of 2 to obtain a tenth result; 

dividing the tenth result by a sensitivity scaling factor; 

recording a positive/negative sign for the tenth result; 

taking the square root of the absolute value of the tenth result to produce a 
eleventh result; 

replacing the positive/negative sign of the tenth result into the eleventh result to 
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produce a twelfth result; and 

incorporating the sixth result and the twelfth result into the seismic profile. 

5. The method of claim 4, further comprising using a second ocean bottom cable 
having a first sensor pair and a second sensor pair, each of the first sensor pair and the 
second sensor pair comprising a pressure sensor and a velocity sensor. 

6. The method of claim 4 or claim 5, wherein each velocity sensor is a geophone. 

7. The method of any one of claims 5 to 6, wherein the pressure sensor is a 
hydrophone. 

8. A method of attenuating reverberations in seismic data, the seismic data including 
pressure data and velocity data containing reverberations, comprising: 

reading the pressure data from at least one pressure sensor; 

reading the velocity data from at least one velocity sensor; 

processing the pressure data and the velocity data to produce composite data 
wherein the reverberations have been reduced, comprising the steps of: 

multiplying the pressure data with the absolute value of the velocity data to 
produce a first result; 

multiplying the velocity data with the absolute value of the pressure data to 
produce a second result; 

summing the first result and the second result to create a third result; 

dividing the third result by a factor of 2 to obtain a fourth result; 

dividing the fourth result by a sensitivity scaling factor; 

recording a positive/negative sign for the fourth result; 

taking the square root of the absolute value of the fourth result to produce a fifth 
result; and 

replacing the positive/negative sign of the fourth result into the fifth result to 
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9. The method of claim 8, wherein the or each velocity sensor is a gebphone 

10. The method of claim 8 or claim 9, wherein the or each pressure sensor is a 
hydrophone. 
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